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We study the asymmetric Hubbard model at half-filling as a generic model to describe the physics 
of two species of repulsively interacting fermionic cold atoms in optical lattices. We use Dynamical 
Mean Field Theory to obtain the paramagnetic phase diagram of the model as function of tem- 
perature, interaction strength and hopping asymmetry. A Mott transition with a region of two 
coexistent solutions is found for all nonzero values of the hopping asymmetry. At low temperatures 
the metallic phase is a heavy Fermi-liquid, qualitatively analogous to the Fermi liquid state of the 
symmetric Hubbard model. Above a coherence temperature, an orbital-selective crossover takes 
place, wherein one fermionic species effectively localizes, and the resulting bad metallic state resem- 
bles the non-Fermi liquid state of the Falicov-Kimball model. We compute observables relevant to 
cold atom systems such as the double occupation, the specific heat and entropy and characterize 
their behavior in the different phases. 
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Quantum degenerate fermions can now be loaded onto 
optical lattices to recreate model hamiltonians with easily 
tunable physical parameters [l[. Two fermionic atomic 
species are minimally required to create a spin-1/2 model, 
where each type of atom is associated to a different pro- 
jection of the spin. This can be achieved by either loading 
the optical lattice with atoms with two different hyper- 
fine spin states or by two different atom species [lr|5|- 
Remarkable aspects of the cold atom systems include the 
variety of lattice-types that can be realized, and the large 
range of interaction strengths that can be accessed, that 
may be up to several times the fermion band width [6|]. A 
lot of recent work has focused on the fermionic Hubbard 
model in 3D optical lattices, which is the paradigm of 
a strongly correlated fermionic system. For instance, a 
BCS-BEC crossover was observed in the attractive case 
using 6 Li atoms 0, [1], and also a Mott insulating state 
was reported in the repulsive case at half-filling using 
40 K atoms [1, [lOl. Although the exact realization of the 
Hubbard model may be achieved in cold atom systems, 
the generic situation, involving two different fermionic 
atomic species, will not exactly fulfill the requirement of 
spin symmetry. Thus, the model typically realized in an 
optical lattice would be an Asymmetric Hubbard model 
(AHM), which interpolates between the Hubbard model 
and the spinless Falicov-Kimball model (FKM), with or 
without population imbalance. This model, in the case 
of attractive interactions, has been considered in recent 
theoretical investigations that focused on the possibility 
of superconducting phases HH. The study of the AHM 
with repulsive interactions is therefore relevant and may 
provide needed guidance for current experimental inves- 
tigations. 



In this paper, we focus on the interesting case of the 
half filled AHM with no population imbalance and the 
associated paramagnetic Mott metal-insulator transition. 
For simplicity, we neglect the harmonic potential in our 
calculations, thus our results apply to the bulk of the 
optical lattice, away from the edges of the trap. We 
study the model within Dynamical Mean Field Theory 
(DMFT), a method which has significantly advanced 
our understanding of the Mott transition phenomenon 
[ill [T3 | . Here, we focus on paramagnetic states which 
may be realized in frustrated optical lattices. Our work 
shows that, akin to the Hubbard case, there is a region 
of coexistent solutions defined by two Mott-transitions at 
critical values U c \ and U C 2 [H| . This region exists within 
the entire asymmetric regime, except right at the FKM 
end-point. While this observation may suggest that the 
AHM qualitatively behaves like the ordinary Hubbard 
model, we find that the actual situation is more subtle 
and cannot be described in terms of a simple interpola- 
tion between the two known limiting cases. In fact, in the 
metallic phase at low temperatures, we find a heavy-mass 
Fermi liquid state which is qualitatively similar to that of 
the spin-symmetric Hubbard model However, above 
a coherence temperature T co h, instead of the standard 
incoherent metal seen in the Hubbard case, an orbital- 
selective crossover occurs, where the solution evolves to 
closely resemble the bad metallic state of the FKM. Thus, 
surprisingly, the behavior of the AHM embodies features 
of both, the Hubbard and the FKM. We also compute 
the double occupancy which is experimentally accessible 
0, and the specific heat and entropy which are useful 
for the estimation of the temperature of experimental 
systems [ic| |. 
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The Asymmetric Hubbard model Hamiltonian reads, 



H 



<ij>;cr=t,4- 



E/^(n* t -l/2)(n 4 -l/2) 



(1) 

c, are the fermion creation and annihilation operators 
and the indices and a label the lattice sites and 

the two species of fermions, respectively. Defining r = 
ti/tf, the AHM interpolates between the FKM at r = 
(where one of the fermion species is immobile) and the 
the spin-symmetric Hubbard model at r = 1. These two 
limits were studied extensively within the framework of 
DMFT (cf. see reviews of Refs. [ll| and HI). We adopt a 
bounded semi-circular density of states for both fermionic 
bands with the half-bandwidths D a = 2t a . This choice 
presents a methodological advantage without affecting 
the qualitative physical behavior of the model [Hj . The 
unit of energy, as customary, is set by = D = 1. 
Within DMFT, the lattice model is map onto a quantum 
impurity problem subject to a self-consistency condition 



imp 
11]. 



The self-consistency condition is given by, 



(2) 



where G a is the local Green's function of the fermions 
with spin a and Qqo- are the bare Green's functions of the 
associated quantum impurity problem fll| . The impurity 
problem can be solved using exact diagonalization and 
quantum Monte Carlo techniques, whose implementation 
has been already detailed elsewhere 

The central quantities in the DMFT method are the 
local Green functions, G a (uj n ), where w„ are the Matsub- 
ara frequencies. Here, since t±we should generically 
have ^ G^. Though the Matsubara Green's functions 
need to be analytically continued to the real frequency 
axis to obtain the spectral function (density of states), 
they nonetheless provide clear physical information. For 
instance, the density of states at the Fermi energy, which 
signals metallic or insulating states, is directly given by 
-l/7rIm[G CT K ->■ 0+)]. 

At half-filling and T = it is well known that the FKM 
has a Mott transition at a critical value Uq KM = D 
whereas the HM is found to have two transitions at 
U% M w 2AD and Ug M « 3D. The transition at U^ M 
is related to the closing of the Mott gap, while that at 
stems from the divergence of the carrier effective 
Both metallic and insulating solutions coexist 
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in the region Uf x M < U < U% M . Given that both the 
Hubbard and FKM have metallic and Mott insulating 
phases, we expect a paramagnetic Mott metal-insulator 
transition in the AHM as a function of the repulsive inter- 
action U for all r. Our results for the Green's functions 
at the generic value r = 0.4, for two values of the in- 
teraction U = D and 3D are shown in the top left and 
right panels of FigQ] The nonzero y-axis intercepts of 
both, G^ and G±, in the data of the left top panel indi- 
cate the presence of (different) metallic states in the two 



channels. On the other hand, as shown in the right top 
panel, at a higher value of the interaction U, the density 
of states vanishes, which signals an insulating state with 
the opening of a Mott gap for the up and down fermions. 
Moreover, as the gap is related to the slope of G a (uj n ) at 
low uj n , we see that the gap has a similar magnitude for 
the two species. Though this establishes the existence 
of a Mott transition, we still need to understand how 
the Mott-transitions of the two end-cases (Hubbard and 
FKM) merge as a function of r and U. 

We first analyze the metallic phase. In the standard 
Hubbard model (r=l), at low temperatures, one has a 
metallic Fermi liquid state with quasiparticles which have 
a renormalized mass to* given by the self-energy £ fl3j | . 
In such a Fermi liquid state, within DMFT, at low fre- 
quencies S(wn) ~ (1 — l/Z)iui n , where Z is the quasipar- 
ticle residue [11] . Inputting the low frequency self energy 
into the definition of the Green's function, we obtain 



G(k,w n ) 



1 



Z 



- e k - (1 - ±)iu). 



(3) 



As §5§ shows, the effective mass renormalization factor 
to* al/X, and at the Mott transition 2^0 and to* — >• 
oo. 

In the asymmetric model it is not a priori evident how 
the interactions will renormalize each one of the band 
masses. The results of the lower right panel of Fig. 2 
show the numerical results for Z a . We find that the mass 
renormalization of the two species is different, with the 
lighter one being more strongly suppressed by the effect 
of interactions. However, they both vanish at the same 
critical point. Thus, the ground-state of the system is an 
asymmetric Fermi liquid with two types of heavy quasi- 
particles, that beyond a critical value simultaneously lo- 
calize into a Mott-Hubbard state. 

As in the Hubbard case, the Fermi liquid state seen in 
the AHM exists below a certain coherence temperature 
T C oh and is characterized by a specific heat which is lin- 
ear in temperature. Above T co /j, the thermal disorder is 
too strong for the quasiparticles to survive and the good 
metallic behavior is gradually lost. Our results for T co h, 
obtained from the specific heat (to be discussed in detail 
later) , are plotted in the lower left panel of Fig. [21 We 
find that this temperature is the same for both species 
consistent with the single Kondo temperature of the as- 
sociated impurity problem ll| . As one increases the tem- 
perature above T co h, the solution has a rapid crossover 
towards a new bad metallic state. Surprisingly, we find 
that this intermediate-temperature metallic state is not 
arbitrary, but essentially coincides with the non Fermi 
liquid solution of the FKM obtained at the same value 
of interaction U. The corresponding numerical result is 
shown in the lower right panel of Fig[TJ where we com- 
pare the solution obtained at T > T co h with the FKM 
solution computed at the same value of T and U. We 
may interpret this rather striking result as if r effectively 
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FIG. 1. (Color online) Matsubara Green's functions of the 
AHM for r=0.4 and T=0.01D. The top left panel shows a 
metallic solution for U =D and the right one show an insulat- 
ing solution for U=3D. The lower left panel shows the pinning 
condition of the two metallic solutions (Im[G(r](un — > 0) = 
2/ D a at T — 0) 11]. The lower right panel shows the same 
Green's functions at a higher temperature T=0.1D > T co h- 
The solution essentially coincides with the FKM solution 
(r=0) computed at the same value of U. 
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FIG. 2. (Color online) The top panel shows the T = phase 
diagram as function of r and U with the Mott transition 
boundaries U& and U C 2- The lower left (right) panel shows 
the behavior of T co h {Z a ) as a function of interaction U for 
various values of the asymmetry factor r. The size of the 
symbols are an estimate of the error bars. 



renormalizes down to zero as T ~ T co h- Since in the FKM 
one of the fermionic species is itinerant while the other is 
fully localized (i^=0), the observed behavior in the AHM 
amounts to an orbital selective crossover to localization 
of the down spin fermions. The AHM can therefore be 
considered as a minimal model to realize an orbital selec- 
tive transition, a concept that has recently received a lot 
of attention in the context of strongly correlated fermions 
in condensed matter systems 3, 15j |. 

The T=0 phase diagram of the model as a function 
of asymmetry r and interaction U is shown in the top 




FIG. 3. (Color online) Double occupancy as function of U 
at T=0.025D, for different values or r=ti/tf. From top to 
bottom, the curves corresponds to decreasing values of r. 



panel of FigfJl As discussed before, at low T, the AHM 
qualitatively behaves as the Hubbard model and has two 
critical Mott transition boundaries, U c .i and U C 2, that 
define a region of coexistent solutions [ll[ . The two crit- 
ical lines merge at the FKM end-point. This coexistence 
survives for a small range of temperatures resulting in a 
first order metal-insulator transition line at low temper- 
atures which culminates in a second order critical point 

111 ] . The calculation of this small critical temperature is 
technically demanding and is left for future work. Sim- 
ilarly to the Hubbard case, we see in Figj2] that for all 
r^O, Z — > as [/ — U C 2 leading to a divergence of the 
effective mass at £/ C 2- 

In order to establish a link between the behavior of 
the model and the physical observables accessible in cold 
atom systems on optical lattices, we compute the double 
occupation (n-j-nj,), which can be used as a thermometer 

161 17 1 . In fact, this quantity is also closely related to 
an order parameter for the Mott metal- insulator transi- 
tion 13, 18j]. The results are shown in Fig |3] for various 
values of U and r at finite T=0.025.D. At larger val- 



ues of r there is a very sharp transition from the metal 
to the insulator state. The metal is characterized by 
linear decrease of the double occupation with increasing 
interaction U, while in the insulator, at larger values of 
the interaction, the double occupation remains low and 
weakly [/-dependent. The rapid variation of (n^n^) with 
a sigmoidal shape between the two regimes is also an 
indication of the finite-temperature region of coexistent 
solutions that we discussed before [18( . At smaller values 
of r, when the system is more asymmetric, the double 
occupation rapidly decreases due to the very small co- 
herence scale and the reduced values of the critical U . 
For the lower r values the temperature of the calculation 
is actually above T co h (cf Figj2j), thus we observe a rather 
slow and wide crossover region for the evolution between 
the metal and the insulator. 

Other experimentally relevant quantities we compute 
are the specific heat and the entropy which are required 
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to estimate the temperature of a cold atom system. In 
fact, these systems are in the sub-millikelvin regime and 
since they do not dissipate energy their evolution is con- 
sidered to occur essentially along isentropic lines [n| . To 
understand the evolution of the specific heat and entropy, 
it is useful to start with the known cases of the Hubbard 
and FKM at T=0. In a regular Fermi liquid the entropy 
vanishes as T — > 0. However, there are also two other 
cases to be considered, which are unusual as they have 
unquenched entropy per site S/N=\n2 at T=0 (we set 
the Boltzmann constant fcs=l). One is the paramag- 
netic Mott-insulator, with one localized particle per site, 
which only has the two spin-states as remanent degree of 
freedom. The other case is the bad metallic state of the 
FKM, where the localized particles are disordered and 
have equal probability of occupying or not each of the 
N lattice sites, thus contributing the In 2 to the entropy 
at half-filling. The itinerant particles, in contrast, adjust 
to the potential created by the former and form a band 
that contributes no entropy at T = 0. Consequently, 
starting from the Fermi liquid state of the AHM, we ex- 
pect strong variations of the specific heat for both, the 
transition to the Mott insulator and the crossover to the 
orbital-selective bad metallic state. 

Our results for the specific heat and entropy are plotted 
in FigHl The left side panels of Fig f4] show the evolution 
of C v in the metallic phase for different values of r. Note 
that for all r ^ and low enough T and U, the system 
is a normal Fermi liquid with the specific heat C v = jT 
and S oc T with 7 proportional to the mean effective 
mass of the quasiparticles. As the temperature increases 
and reaches T co h, the singlet Fermi-liquid metallic state 
breaks down and the system crosses over to the orbitally 
selected bad-metal, which has FKM-like behavior, thus 
having an entropy of In 2. This is clearly seen in the lower 
right panel of FigHJ where after an initial linear increase, 
the entropy attains a small plateau around In 2. This 
plateau is concomitant with the first maximum observed 
in the C V (T) at a T ~ T co h- Physically, it corresponds 
to the orbital selection, i.e, the effective localization of 
the heavy-particle band. As T is further increased, the 
remaining degrees of freedom will eventually contribute 
their In 2 to the entropy, producing a second maximum in 
C v at the bare energy scale T ~ t^=D/2. In the high-T 
limit the specific heat decreases and the entropy saturates 
to the maximal value In 4. For comparison, the behavior 
of the specific heat and the entropy of a typical Mott 
insulating state at higher interaction U — 2D and r = 0.4 
is also shown in the right panel of FigfJJ At low T, C V (T) 
shows activated behavior indicating the presence of the 
Mott gap. We also observe, consistent with our previous 
discussion, that the entropy starts from the unquenched 
value In 2, and reaches In 4 in the high-T limit. 

To conclude, we have studied the behavior of the asym- 
metric Hubbard model as a generic model for two species 
of fermionic cold atoms loaded onto an optical lattice. 
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FIG. 4. Left panels: Specific heat as function of T for U—D 
and different values of r. For r=l the linear Fermi-liquid 
behavior is observed. As r decreases, a two-peak structure is 
visible. The first peak moves to lower T with decreasing r, 
revealing the reduction of T co h and the increase of the effective 
mass (from the increase of the slope) .Right panels: C v (top) 
and entropy (bottom) as a function of T for r=0A and U=D 
(metal) and U=2D (Mott insulator). 

We determined the phase diagram and found a Mott 
metal-insulator transition for all degrees of asymmetry. 
The Mott transition scenario is the same as that of the 
Hubbard model for all nonzero values of the asymme- 
try r. Within the metallic state, we found and inter- 
esting temperature-driven orbital selective crossover. At 
low temperatures the AHM has the same qualitative be- 
havior as the symmetric Hubbard model, however with 
two different heavy mass Fermi-liquids. These two liq- 
uids have, nevertheless, a single coherence temperature. 
Above a coherence temperature, one fermionic species lo- 
calizes and the system exhibits non-Fermi liquid behav- 
ior similar to that of the Falicov-Kimball model. We also 
computed the double occupation and some thermody- 
namical observables which may permit the experimental 
identification of this remarkable physical behavior. 
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